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ABSTRACT 

We report on the first Green Bank Telescope observations at 427, 820 and 
1400 MHz of the newly discovered, highly inclined and relativistic double pulsar 
binary. We focus on the brief eclipse of PSR J0737— 3039A, the faster pulsar, 
when it passes behind PSR J0737— 3039B. We measure a frequency-averaged 
eclipse duration of 26.6 ± 0.6 s, or 0.00301 ± 0.00008 in orbital phase. The 
eclipse duration is found to be significantly dependent on radio frequency, with 
eclipses longer at lower frequencies. Specifically, eclipse duration is well fit by a 
linear function having slope (—4.52 ±0.03) x 10~^ orbits MHz~^. We also detect 
significant asymmetry in the eclipse. Eclipse ingress takes 3.51 ± 0.99 times 
longer than egress, independent of radio frequency. Additionally, the eclipse 
lasts (40 ± 7) X 10~^ in orbital phase longer after conjunction, also independent 
of frequency. We detect significant emission from the pulsar on short time scales 
during echpse in some orbits. We discuss these results in the context of a model 
in which the eclipsing material is a shock-heated plasma layer within the slower 
PSR J0737— 3039B's light cylinder, where the relativistic pressure of the faster 
pulsar's wind confines the magnetosphere of the slower pulsar. 
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1. Introduction 

The discovery of the double-pulsar binary system PSR J0737— 3039 (Burgay et al. 2003; 
Lyne et al. 2004) represents a landmark in neutron-star astrophysics. This relativistic system 
consists of two neutron stars, both radio pulsars, in a 2.4-hr orbit having eccentricity 0.08, 
viewed nearly edge-on. The binary promises unparalleled tests of General Relativity. During 
each orbit, the faster- rotating of the two pulsars, PSR J0737— 3039A (hereafter A), which 
has period P — 22 ms, is eclipsed briefly at conjunction, while the slower-rotating of the two, 
PSR J0737— 3039B {P — 2.7 s; hereafter B), has unusual and longer- lived flux enhancements 
near the same epoch (Lyne et al. 2004). This bizarre and unprecedented behavior offers a 
new opportunity to understand pulsar magnetospheric and wind structures. 

This paper is part of a short series describing the first follow-up observations of this 
remarkable source using the 100-m diameter National Radio Astronomy Observatory Green 
Bank Telescope (GBT).^ Here we concentrate on the properties of the pulsar "A" eclipse. 
From observations at the Parkes and Lovell Telescopes, Lyne et al. (2004) have already 
established that the eclipse has duration 20-30 s and is roughly independent of frequency. 
The greater sensitivity of the observations described here reveal that the eclipse duration is 
in fact slightly but significantly dependent on radio frequency, as well as demonstrate for the 
first time that the eclipse is asymmetric. 

2. Observations 

Our observations with the GBT were promptly scheduled as a NRAO Rapid Response 
Science Program^ four times in 2003 December and once in 2004 January and made use of 
receivers at 427, 820 and 1400 MHz. The observational details are summarized in Table 1. 
The 427 and 820 MHz receivers were located at the prime focus of the GBT, are cooled 
FET amplifiers, and have approximate system temperatures of 57 K and 25 K, respectively. 
The 1400 MHz receiver was located at the Gregorian focus and is a cooled HFET amplifier 
with system temperature of 20 K. All provide one beam and are sensitive to orthogonal 
polarizations. Signals from the receivers are transferred from the telescope to the control 
room via optical fibre. 

The data reported on here were acquired with a single Berkeley-Galtech Pulsar Machine 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 

^ http : / / www . via. nr ao .edu/astro/prop/rapid/ 
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(BCPM). The BCPM is an analog/digital filter bank which samples 2 x 96 channels using 
4 bits at flexible sampling rates and channel bandwidths (Backer et al. 1997). For our 
observations, the two orthogonal polarizations were summed in hardware prior to being 
written to disk. The sampling time for all observations was 72 /is. The total bandwidth 
used (and hence the channel bandwidth) depended on the observing frequency (Table 1). 

3. Analysis and Results 

The data were analyzed using the PRESTO software package (Ransom 2001). We dedis- 
persed the data using dispersion measure 48.9115 pc cm~^ and did the following analysis on 
the resulting time scries. First, we folded data obtained within 2 min of pulsar A's eclipse 
with the cphcmcris reported by Lyne et al. (2004) using 64 phase bins. The resulting profile 
was cross-correlated with a high signal-to-noise template (obtained using data from several 
complete orbits at the appropriate radio frequency) and the temporal offset was recorded. 
Then, we folded data from within the 2 min of the echpse again, but in 2-s intervals. We 
then determined the pulsed intensity in each 2-s interval, first by aligning the pulse with 
the template using the previously recorded offset, then by finding the best-fit scale factor 
between the template and aligned data profile. The best-fit value was found by mini- 
mization, and the equivalent la uncertainties determined from the values of the scale factor 
corresponding to the minimum plus 1. Note that no effort was made to calibrate the in- 
tensities in physical units; we assume that the telescope gain and system temperatures were 
constant throughout each observation. The average pulse arrival time in each 2-s interval 
was recorded and reduced to the solar system barycenter using the published position for 
the pulsar (Lyne et al. 2004). In this way, we obtained a pulsed flux intensity time series 
consisting of a bary centric time and an intensity for each of the 120 2-s intervals around 
conjunction. The number of conjunctions observed at each frequency is given in Table 2. 
Note that we found no evidence for any changes in dispersion measure near eclipse. Us- 
ing the 427 MHz data, we set an upper limit on any dispersion- measure variation on > 2-s 
time scales of 0.016 pc cm~'^. This was done by running the pulsar timing software package 
TEMPO^ on the 427 MHz data (with eclipse points identifled by large arrival-time uncertainty 
and omitted) without fltting for any parameters, and taking the largest residual as an upper 
limit on any delay due to an anomalous dispersion-measure change. 

Average eclipse light curves for each frequency were obtained by combining results of 
every observed conjunction and arc shown in Figure 1. The light curve at 427 MHz has more 
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scatter than the other hght curves, primarily due to higher levels of sporadic, short-time scale 
radio interference of unknown origin, particularly during the first orbit. 

We detect significant variations in the pulsed emission at the same orbital phase but in 
different orbits, during the eclipse. An example of such variability is illustrated in Figure 2, 
in which the light curves for the three 1400 MHz echpses are shown. In the second observed 
orbit, there is an apparent 5(t detection of the pulsar during the nominal eclipse period, as 
indicated in the Figure. Investigation of this 2-s interval's folded profile by eye confirmed 
that the pulse, with the expected pulse shape and phase (though appropriately noisy), was 
present (see Fig. 2, inset). Also, in Figure 1, there is emission in the eclipse region both before 
and after conjunction in the average 427 MHz light curve. Examination by eye confirmed 
the approximate pulse shape at the expected phase in several individual 2-s intervals in two 
of the three 427 MHz orbits. No such pulse appearance was seen in any of the 820 MHz 
eclipses, however. 

Several new properties of the eclipse are apparent from Figure 1. First, the eclipse is 
clearly asymmetric in that the eclipse ingress is longer in duration than egress. Also, the 
eclipse appears to last longer post conjunction. 

To quantify these properties, we fit the orbit- averaged light curves at each frequency 
with two functions of the form 

Here -F(0) is the flux at orbital phase 0, defined such that conjunction is at = 0. This 
function, having two free parameters 0o and w, has value unity for — (f)o)/w << — 1, zero 
for {(f) — (f)o)/w » 1, and one half for (f) — (pQ. The transition from values near unity to 
zero occurs in a width ~ 2w. Optimal values of 0o a-nd w for eclipse ingress and egress were 
determined separately for each light curve in Figure 1 using straightforward minimization 
and data in the interval —0.006 < < for ingress and < < 0.006 for egress. With two 
free parameters per side, the number of degrees of freedom per fit, u, was 23. The best fits 
are shown as solid curves in Figure 1; they generally yield values of reduced ii^ar unity, 
except for the 427 MHz light curve, for which the higher results from the extra scatter 
noted above. Contours of constant x^ <Po/w space were approximately circular, indicating 
that these parameters were largely independent. Best- fit parameters and xl values are given 
in Table 2. 

Figure 3 summarizes the results of these fits. Panel (a) shows the echpse FWHM, 
defined as the sum of ingress and egress values of 100 1. The clear trend in the data indicates 
that the echpse duration is frequency-dependent. Indeed the reduced of a fit to the 
data assuming a simple constant mean value is 6.1, clearly ruling this model out. The data 
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are very well fit by a linear model, having slope (—4.52 ± 0.03) x 10~^ orbits MHz~^, and 
intercept (3.412±0.003) x 10~^. This predicts that the eclipse will disappear at approximately 
7.5 GHz, assuming no other eclipse mechanism becomes important. The mean eclipse FWHM 
duration, averaged over these three frequencies, is 0.00301 ± 0.00008 in orbital phase, or 
26.6 ±0.7 s. 

Figure 3, panel (b), shows ingress and egress values of w. One aspect of the asymmetry 
of the eclipse is very apparent from this plot: all ingress values of w are higher than the egress 
values. Thus, the pulsar exits eclipse faster than it enters. We quantify this effect by taking 
the ratio of the ingress to egress values of w, Wi/we- We find, for 427, 820 and 1400 MHz, 
Wi/we = 3.34 it 0.59,3.74 ± 0.78, and 3.45 ± 1.40, respectively. Thus, this aspect of the 
eclipse asymmetry is not dependent on frequency. The mean Wi/we is 3.51; the uncertainty 
on this value, estimated from the uncertainties at each frequency, is 0.99. Estimating the 
uncertainty from the RMS scatter in the three values suggests a much smaller uncertainty 
of 0.15; this may be a chance occurence due to small- number statistics. 

Figure 3, panel (c), shows the fitted values of 0o for ingress and egress at each fre- 
quency. Clearly, |0o| values for egress are systematically larger than for ingress (Fig. 2c), 
demonstrating that the eclipse lasts longer post-conjunction. Note that the latter is opposite 
to what would be expected if the w asymmetry were the only one present, and hence is a 
different effect. We quantify this using the difference between ingress and egress values of 
00- For 427, 820 and 1400 MHz, we find differences of 0.00042 ± 0.00007, 0.00036 ± 0.00006, 
and 0.00042 ± 0.00009, respectively. Thus we find no evidence for frequency dependence of 
this eclipse asymmetry. That the absolute values of 0o for both ingress and egress decrease 
significantly with increasing frequency is a restatement of the fact that the eclipse duration 
decreases with increasing frequency. 

4. Discussion 

The best-fit inclination angle of the system as reported by Ransom et al. (2004) is 
88.7° ± 0.9°. Given the measured projected semi-major axis of pulsar A's orbit, this implies 
that A's beam must pass within 0.07 It-s of pulsar B. The latter has light-cyhnder radius 
Ric — cP/27r = 0.45 It-s. Thus, if B's magnetosphere were unperturbed, then A's beam as ob- 
served by us would have passed well within its light cylinder and through its magnetosphere. 
However, as discussed by Lyne et al. (2004), given the relative spin-down luminosities of the 
two pulsars, with A's being a factor of ~3625 greater than B's, it seems inevitable that A's 
wind has disrupted B's outer magnetosphere at least halfway inside B's light cylinder and 
possibly further out than A's beam passes at conjunction. At A's conjunction, the relative 
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transverse velocities of the pulsars is 680 km s"^, so the eclipse duration implies an approxi- 
mately constant physical size of the eclipsing region of 18,100 km, or 0.060 It-s, comparable 
to the estimated 0.07 It-s impact parameter, though in a direction parallel to the orbital 
plane. This clearly implies an eclipse region that is much smaller than B's unperturbed 
magnetosphere. 

Synchrotron absorption in a shock heated "magnetosheath" surrounding and contain- 
ing B's magnetosphere is a possible model for the A eclipse (Arons et al. 2004a; Lyutikov 
2004, Arons et al., in preparation). A's wind probably confines B's magnetosphere on the 
side facing A (which changes periodically with B's rotation over 2.7 s) such that the latter 
likely resembles a time-dependent variant of the Earth's magnetosphere: compressed on the 
side facing the source of the wind (the Sun in the terrestrial case, A in this case) with a 
magnetotail on the opposite side. The plasma from A passes through a bow shock which 
spreads and weakens with distance down the flanks of the tail. The wind plasma from A 
which passes through the bow shock, flows in a magnetosheath layer between the bow shock 
and the magnetopause current layer that confines B's magnetosphere. If the bow shock 
compresses and heats the A pulsar's wind plasma as if it occurred in a weakly magnetized, 
semi-infinite fiow, then synchrotron absorption in this shock heated plasma can account for 
the strong eclipse of A by B, albeit for surprisingly high pair densities in A's wind (Arons 
et al. 2004a). The dimensions of the magnetopause and bow shock on the side facing A, 
along with the impact parameter of the line-of-sight with respect to B at conjunction, are 
such that A's radio beam can pass through the northernmost or southernmost extremities 
of the magnetosheath layer (where north or south means with respect to the orbital plane.) 

Such a scenario is consistent with the sharp edge and mildly frequency-dependent shape 
of the eclipse, since the synchrotron absorption optical depths at the observation frequencies 
used in this study are large and the shock causes the plasma density and relativistic temper- 
ature to jump up sharply across the thin shock layer. The eclipse-to-eclipse variations and 
occasional appearances of the pulse during the nominal eclipse region suggest a changing 
eclipse medium. If pulsar B is an oblique rotator, it is possible that the rotation of its mag- 
netosphere causes the magnetospheric and bow shock shapes to be asymmetric on average 
between ingress and egress by amounts which account for the observed asymmetries. 

If this picture and physical model is correct, one reaches the interesting conclusion that 
A's wind shocks as if it is not strongly magnetized at a distance of only ~750 times its 
light cylinder radius. Previous conclusions regarding pulsar wind magnetization have been 
derived from the behavior of relativistic winds as they pass through their termination shocks 
in pulsar wind nebulae (see Kaspi et al. 2004, for a recent review), located at much larger 
distances (typically ~ 10^ fight cyfinder radii) from their driving pulsars (e.g. Kennel & 
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Coroniti 1984). 

The above A eclipse model (see Arons et al. , in preparation) predicts that at sufficiently 
high radio frequencies, the eclipse will clear, consistent with the observed frequency depen- 
dence of the eclipse duration. For nominal, highly simplified models, transparency might 
appear at frequencies above 5-10 GHz, roughly quantitatively consistent with the observed 
frequency dependence extrapolation. If the impact parameter of the line-of-sight with re- 
spect to B during eclipse is sufficiently small, cyclotron scattering within B's magnetosphere 
might also contribute to the high-frequency opacity. However, nominal values for the geom- 
etry suggest this to be unlikely. In addition, the eclipse duration and morphology should 
change with time on few-year time scales as relativistic periastron precession, and, perhaps 
more importantly, geodetic precession, which is expected to be very important in this system 
(Lyne et al. 2004), change the alignment of the pulsars. This will alter both our viewing 
vantage point and the dynamic pressure of A's wind on B at conjunction. 

The durations of pulsar A ingress and egress are comparable to the rotational period 
of B. As B's magnetosphere could well present a different response to A's wind pressure 
depending on B's rotational phase, significantly different A eclipse behavior from orbit to 
orbit might be expected. A higher time resolution investigation with higher sensitivity may 
reveal B phase dependent effects on the eclipse. 

We pay tribute to the Parkes survey team for so wonderful a discovery. We thank C. 
Bignell, F. Ghigo, and G. Langston for assistance with the GBT observations, and M. Roberts 
and B. Rutledge for useful discussions. This work was funded by an NSERC Discovery Grant 
and Steacie Fellowship Supplement to V.M.K. Additional funding to V.M.K. was provided 
by FQRST and CIAR. J.A. was partially supported by NASA grants TM4-5005X, NAG5- 
12031, G03-4063B and HST-HF-01157.01-A. He also thanks the taxpayers of Cahfornia for 
their continued indulgence. A.S. acknowledges support provided by NASA through Chandra 
Fellowship PF2-30025 awarded by the Chandra X-ray Center, which is operated by the 
Smithsonian Astrophysical Observatory for NASA under contract NAS8-39073. 

Note added in proof. - Very recently, McLaughlin et al. (2004) report modulation of 
A's flux at B's period during eclipse ingress, from a higher time resolution analysis of the 820 
MHz data set used here. This suggests that the eclipse medium is indeed highly dependent 
on B's rotational phase. 
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Table 1: Summary of GBT Observations of PSR J0737-3039 



Date 


Start UTC 


Start MJD 


Duration 


Frequency 


Bandwidth 


No. Channels 


((l(l/iiiiii/yy) 






(hr) 


(MHz) 


(MHz) 




11/12/03 


04:36:27 


52984.19199 


6.15 


1400 


96 


96 


18/12/03 


04:12:15 


52992.17519 


5.30 


820 


48 


96 


23/12/03 


05:22:22 


52997.22388 


4.72 


820 


48 


96 


01/01/04 


03:30:09 


53005.14595 


6.22 


427 


48 


96 



Table 2: Fit Parameters for the PSR J0737-3039A EcUpse 



Frequency 


No. Eclipses 




Ingress 






Egress 






Observed 


00 


w 


X23 


00 


w 


X23 


(MHz) 




(xlO-^) 


(xlO-3) 




(xlO-3) 


(xlO-3) 




427 
820 


3 
4 


-1 60+°°^ 

_1 Q4+0.04 


"Sfi+o oe 

>J.OD_Q05 

n 40+0.05 

^•^^-0.04 


2.18 

0.834 


1 or,+0.03 
^■°^-0.08 

-'-•''"'-0.03 


A -1 7+0.03 
'"'■-'-'-0.02 

1 3+°°3 
'"'■-'-'-'-0.02 


1.65 
0.76 


1400 


3 


1 -1O+0.08 
-i-i»_0.06 


ci-^+o os 
u.oo_o.o8 


0.53 


1 60+°°^ 
^■'-''J-0.04 


'^■^"J-0.06 


1.02 



Note: X23 is the reduced for 23 degrees of freedom. 
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Orbital Phase (w.r.t. superior conjunction) 



Fig. 1. — Pulsar A eclipse light curves. Each point represents 2 s of data; the shown curve 
is 4-min in duration, centered on conjunction. The x-axes are orbital phase with respect to 
conjunction. The y-axes are pulsed flux, normalized such that the pre-eclipse flux is unity. 
The panels are for (a) 427 MHz, (b) 820 MHz, and (c) 1400 MHz. The solid vertical line 
indicates conjunction and the horizontal dotted lines show flux. Vertical dashed lines at 
= —0.006, +0.006 indicate the range of data fltted. The best-flt ingress and egress model 
curves are shown as solid lines. Best-flt parameters are provided in Table 2. 
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Fig. 2. — Close-up of eclipse region for the three orbits observed at 1400 MHz. The first, 
second and third orbits are indicated with circles, triangles and stars, respectively. Note the 
circled second orbit point just before conjunction; the pulsed intensity of this 2-s integration 
is 5cr above zero, indicating significant pulsed emission during eclipse. The pulse profile for 
this 2-s interval is shown in the inset (solid line), as is the average template profile (dotted 
line), with the pre-determined phase alignment. 
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Fig. 3. — Various fitted properties of pulsar A's eclipse as a function of radio frequency. Panel 
(a): eclipse duration (defined as the sum of the ingress and egress values of |0o|, equivalent 
to FWHM), and best-fit linear model (solid line). Panel (b): values of w (approximately half 
the time for transition into and out of eclipse) for ingress (open symbols) and egress (filled 
symbols). Panel (c): values of \4>o\ for ingress (open symbols) and egress (filled symbols). 



